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ABSTRACT 

We employ the A-variance analysis and study the turbulent gas dynamics of sim¬ 
ulated molecular clouds (MCs). Our models account for a simplified treatment of 
time-dependent chemistry and the non-isothermal nature of the gas. We investigate 
simulations using three different initial mean number densities of no = 30,100 and 
300 cm“^ that span the range of values typical for MCs in the solar neighbourhood. 
Furthermore, we model the CO line emission in a post-processing step using a radiative 
transfer code. We evaluate A-variance spectra for centroid velocity (CV) maps as well 
as for integrated intensity and column density maps for various chemical components: 
the total, H 2 and ^^CO number density and the integrated intensity of both the ^^CO 
and ^^CO (J = 1 ^ 0) lines. The spectral slopes of the A-variance computed on the 
CV maps for the total and H 2 number density are significantly steeper compared to 
the different CO tracers. We find slopes for the linewidth-size relation ranging from 
0.4 to 0.7 for the total and H 2 density models, while the slopes for the various CO 
tracers range from 0.2 to 0.4 and underestimate the values for the total and H 2 density 
by a factor of 1.5 — 3.0. We demonstrate that optical depth effects can significantly 
alter the A-variance spectra. Furthermore, we report a critical density threshold of 
^ 100 cm“^ at which the A-variance slopes of the various CO tracers change sign. We 
thus conclude that carbon monoxide traces the total cloud structure well only if the 
average cloud density lies above this limit. 
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1 INTRODUCTION 


The interstellar medium (ISM) is dominated by highly tur¬ 
bulent motions, which contribute to regulating stellar birth 


in molecular clouds (MCs) (see, e.g. Mac Low & Klessen 
2004| [Scalp V Elmegreen||2004| |Elmegreen V Scalo 2004| 


McKee V Ostriker|2007| Hennebelle V Falgarone|2012 ). On 

large scales, turbulent gas motions support MCs against 
gravitational collapse. However, the shocks associated with 
supersonic turbulence will create overdense regions on small 
scales, which in turn may collapse and form stars and clus¬ 
ters. In addition, ISM turbulence influences other physical 
processes such as the chemical makeup of the gas and the 
associated heating and cooling processes, the efficiency with 
which the external radiation field will be able to penetrate 
into the interior of dense clouds, as well as the overall mag¬ 
netic field. Better understanding the role of turbulence in 
the ISM, therefore, is of pivotal importance for many fields 
of modern astrophysics (see also the lecture notes by Klessen 
fc Glover|20l4| . 


Comparing theoretical simulations to observational 
measurements is very difficult, since observations are always 
a complex convolution of the density and the velocity field, 
affected by several other important physical processes (e.g. 
by magnetic fields, the chemistry, stellar feedback, etc.). 
Furthermore, it is often difficult to infer reliable physical 
parameters from observational measurements. For example, 
observers rely on abundant tracers, e.g. ^^CO, in order to 
measure the amount of H 2 gas in MCs. An obstacle with the 
line emission of ^^CO is that it becomes optically thick in 
dense cloud regions and thus is not a good tracer for those 
regions anymore. To avoid this problem, observers also em¬ 
ploy other tracers to study the mass distribution within a 
cloud, e.g. dust or ^^CO, which is often optically thin. In 
addition, molecular tracers like CS, HCN, HCO^ or NH 3 
can be used to study the gas mass in very dense cloud re¬ 
gions. However, comparisons with simulations are needed in 
order to study the influence of chemical inhomogeneities and 
optical depth effects on observational measurements. 


Lis et al. (1996) introduced maps of centroid velocities 
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(CV) and showed that those are to a certain degree sen¬ 
sitive to the underlying physics. Accordingly, it is possible 
to use two-point statistics in order to recover important in¬ 
formation (e.g. optical depth effects, density and velocity 
fluctuations, etc.) of an astrophysical system, e.g. by us¬ 
ing centroid velocity increment structure functions (see, e.g. 
Hily-Blant et al.|2008 Federrath et al.|2010 ), spectral corre¬ 
lation functions ( Rosolowsky et al.|1999 ), the velocity chan- 


2 METHODS AND SIMULATIONS 

The simulation data presented in the following sections are 
also used in Bertram et al.|(|2014| ) and Bertram et cir| (2015a). 
While Bertram et al. (2014) use the technique of Principal 


nel analysis (Lazarian & Pogosyan 2000 2004), the Prin¬ 


cipal Component Analysis (Heyer & Schloerb 1997[ Brunt 


& Heyer 2002a|) or the A-var iance (see, e.g. Stutzki et al. 


1998| |Mac Low fc Ossenkopf||2000| |Ossenkopf et al.||200l[ 


Benscti et al.||2001| [Ossenkopf et al.||2008a|b ). The latter 

provides a wavelet-based measure for characterizing struc¬ 
tures in astronomical datasets. For example, [Stutzki et al.| 
(1998) used the A-variance in order to measure the scaling 


Component Analysis (PCA) in order to statistically ana¬ 
lyze the turbulent flows in the spectral data cubes, [Bertram] 
et al. [ (2015a) evaluate the slopes of the structure functions of 
2D projected centroid velocities and compute Fourier spec¬ 
tra. Here, we summarize the most important aspects of our 
hydrodynamical simulations and the radiative transfer post¬ 
processing. 


2.1 Computational method 

The simulations in this paper were performed using a mod- 


behavior of structures in observed ima ges by analyzing the ifigd version of the Zeus-Mp MHD code ||Norman| [2000 
line emission of and ^®CO, while 


Bensch et al. 


( 2001 ) 


studied the influence of white noise and beam smoothing on 


the A-variance spectra. Moreover, Ossenkopf et al. (2008b) 


applied the A-variance to interstellar turbulence and obser¬ 
vations of MCs and tested the capabilities of the method 
in a practical use by applying different filter functions with 
different diameter ratios. 


Several studies tried to reveal the influence of pro¬ 
jection effects on different statistical quantities, finding 
that various physical processes might influence the projec¬ 
tion of three-dimensional data on a two-dimensional map 
on the sky. For example 


Lazarian Pogosyan (2004), 
Burkhart et al.[ (2013b), Burkhart et al. (2014) and Bertram 


et al.[ (2015a) have shown that optical depth effects might 


significantly alter the statistics of centroid velocities due to 
different opacities of the gas tracers. Furthermore, the tur¬ 
bulent driving as well as density and temperature fluctua¬ 
tions along the line-of-sight also have a significant impact 
on the CV statistics ( Lazarian V Esquivel[[2003 Ossenkopf 


et al. 2006 Esquivel et al. 2007| Hily-Blant et al. 2008). 


Much progress has been made during the last years in this 
context, although we are still missing a consistent picture of 
turbulence theory in the ISM. 

In this paper we study high-resolution, 3D and time- 
dependent chemistry models of hydrodynamical numerical 
simulations of the turbulent ISM. We use the A-variance 
method (see, e.g, 


Stutzki et al. 1998 Bensch et al. 2001 


Ossenkopf et al. 2008a) in order to study the structure of 


MCs and analyze how different chemical tracers affect the 
A-variance spectra. We also vary the initial number den¬ 
sity in order to explore how different densities affect the 
statistics. Moreover, we perform a radiative transfer post¬ 
processing and produce synthetic maps as well as position- 
position-velocity (PPV) cubes of the ^^CO and ^^CO emis¬ 
sion. We then compute maps of centroid velocities, column 
densities and integrated intensities in order to analyze how 
the chemical inhomogeneity and the variable opacity affect 
the various A-variance spectra. 

In Section!^ we present our numerical simulations, in¬ 
troduce the radiative transfer post-processing and explain 
the statistical methods. In Sectionwe present our results, 
which are discussed in Section Finally, we summarize our 
findings and present our conclusions in Section 


Hayes et al. 2006). We have embedded a detailed atomic 
and molecular cooling function, described in [Glover et al.[ 
(2010) and Glover & Glark (2012), together with a simpli¬ 


fied treatment of the molecular gas chemistry. The chemical 


network is based on the work of Nelson V Langer (1999) 


and Glover V Mac Low (2007), and allows us to follow the 


formation and destruction of H 2 and GO self-consistently 
within our simulations. The network tracks the abundances 
of 9 species and follows 30 chemical reactions. We adopt 
the standard solar abundances of hydrogen and helium. The 
abundances of carbon and oxygen are taken from [Sembach] 


et al. 


(|2000l. We use xc = 1-4 x 10“^ and xo = 3.2 x 10“^ 


where xq and xq are the fractional abundances by number 
of carbon and oxygen relative to hydrogen (for further read¬ 
ing see also Glover et al. 2010 Glover V Mac Low 2011). 


When we start our simulations, the carbon is assumed to 
be singly ionised as G^, while the hydrogen, helium as well 
as the oxygen are in atomic form. For further information 
about the chemical model we refer the reader to I Glover ^ 


Glark (2012). 


We run simulations with three different initial num¬ 
ber densities no of the hydrogen nuclei, no = 300,100 and 
30 cm“^. The temperature of the gas at the beginning of 
our runs is set to a constant value of 60 K. We use an uni¬ 
form weak magnetic field strength Bq = 5.85 //G, which is 
initially oriented parallel to the z-axis of the computational 
domain. We do not include self-gravity. The gas is initially 
uniform and embedded in a periodic box with a side length 
of 20 pc. The turbulent simulations are uniformly driven be- 
tween wavenum bers k — 1 and k — 2 ( Mac Low et al.[|1998 


Mac Low]|l999 ) with a 3D rms velocity of r'rms = 5 km s ^. 


This value remains approximately constant throughout the 
whole simulation period. For the dust-to-gas ratio, we adopt 
the standard local value of 1:100 ( Glover et al.|201(3 ), and as¬ 
sume that the dust properties do not vary with the gas den¬ 
sity. The cosmic ray ionization rate was set to C = 10“ 

For the incident ultraviolet radiation field, we adopt the 
standard parameterization of Draine (1978). This field has 


a strength Gq — 1.7 in Habing (1968) units, corresponding 
to an integrated flux of 2.7 x 10“^ ergcm“^s“^. We use a 
numerical resolution of 512^ grid cells and discuss the in¬ 
fluence of our limited numerical resolution in more detail in 
Appendix [A| Using this value for the resolution as well as 
the scale of the total box, we can compute the cubic cell size 
and obtain Ax ~ 0.04 pc for our numerical simulations. 
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2.2 Radiative transfer post-processing 

We use the radiative transfer code Radmc- SeQ ( [Dullemond 
2012 ) in order to model the CO (J = 1 ^ 0) line for both 


and ^^CO. Furthermore, we apply the Large Velocity 
Gradient (LVG) approximation ( Sobolev]|1957|) to compute 
the level populations, as explained by Shetty et ah ( 2011 a). 
Beyond the line emission, Radmc-3d also accounts for the 
dust continuum emission, which we subtract off before we 
analyze the statistical properties of the synthetic maps. 
We use a number of 512 channels in velocity space for 
our radiative transfer post-processing, corresponding to a 
spectral resolution of ~ 0.07kms“^, ~ 0.06 kms“^ and 
~ 0.05 kms“^ for the n300, nlOO and n30 model, respec¬ 
tively. 

Our simulations do not explicitly track the abundance 
of (which would be costly), and so we need a proce¬ 

dure to relate the ^^CO number density to that of ^^CO. 
A common assumption is that the ratio of ^^CO to ^^CO is 
identical to the elemental abundance ratio of to (see, 
e.g. Roman-Duval et al.|201Q ). In the majority of our analy¬ 
sis, we make the same assumption and set the ^^CO to ^^CO 
ratio to a constant value, R 12/13 = 50, which we use to gen¬ 
erate a first set of GO column density maps. However, the 
effects of chemical fractionation ( [Watson et al.|~976|) and 


selective photodissociation of ^^GO (see, e.g. 


Visser et al. 


2009) can significantly alter the value of R 12/13 wi thin the 
cloud ( Rollig V Ossenko^|2013 Sziics et 'ar]|2014 ). There¬ 
fore, we also explore the effect that this may have on the 
statistics of the A-variance, using numerical results of jSzucsj 
et al. (2014). Sziics et al. (2014) give a numerical htting for¬ 
mula, which relates the ratio R 12/13 to the ^^GO column 
density in the cloud. Thus, we also produce a second set of 
maps using their routine and compute variable ^^GO num¬ 
ber densities by dividing the ^^GO number densities by the 
mean value of R 12/13 from the fitting formula for each line- 
of-sight. Finally, we can use these ^^GO number densities 
in order to compute the ^^GO emission in the same fash¬ 
ion as in our constant ratio models. The differences in the 
A-variance spectra of the two sets of ^^GO emission maps - 
those derived using a constant R 12/13 and those that use a 
spatially varying value of R 12/13 - will be discussed in Sec¬ 
tion |3]3] below. 

The radiative transfer calculation yields position- 
position-velocity (PPV) cubes of brightness temperatures 
Tb, which are related to the intensity via the Rayleigh-Jeans 
approximation. 


r,M = (£) 


VJ 2fcB’ 


( 1 ) 


where is the specific intensity at frequency v and ks the 
Boltzmann constant. We will refer to these as the “inten¬ 
sity” models. In analogy, we construct centroid velocity and 
column density maps out of the PPP simulation data of 
the density and velocity held, as described in the following 
Section [23] We will refer to these as the “density” models, 
following the notation already used in Bertram et al. (2014) 


and Bertram et al. (2015a). 


www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d/ 


L[pc] 



L/D 

Figure 1. Example of a A-variance spectrum, plotted as a func¬ 
tion of spatial scale, normalized by the total box size and averaged 
over all available time snapshots. The inferred slope cx is also in¬ 
dicated (solid line). The A-variance is computed on a CV map, 
where the velocities are weighted by the total density held for an 
initial number density of no = 300 cm“^. The htting range is in¬ 
dicated by a horizontal dashed line. Error bars denote 1-cr spatial 
and temporal variations. 


2.3 Centroid velocity, integrated intensity and 
column density maps 

One part of our study is based on the statistics of centroid 
velocities. The centroids C(x, y) — C(y) are dehned as 


C(r) 


f F(r,z)vz(r,z)dz 
f F(r, z)dz 


( 2 ) 


where the variable Vz(r,z) is the velocity component along 
the line-of-sight (e.g. the ^-direction). Thus, C(r) is a map 
of line-of-sight projected velocities. The quantity F(r,z) is 
a statistical weight. It denotes either the underlying density 
held or the brightness temperatures from the PPV cubes. As 
shown by Lis et al. (1996), centroid velocity maps are very 
sensitive to the underlying physics of a MG, which means 
that they can be used in order to infer important information 
about the astrophysical system. Furthermore, GV maps can 
be easily computed from observational data that rely on 
spectral measurements. These can be then translated into 
maps of centroid velocities via equation §. 

Beside the GV maps, we also consider maps of inte¬ 
grated intensity and column density. We calculate the for¬ 
mer along a given line-of-sight (e.g. the ^-direction) via 



TB(r,z) dv, 


(3) 


where Tb is the brightness temperature, as introduced in 
Section [ 2 ?^ The latter is given by 



n(r, z) dz, 


(4) 


where n is the number density of an individual chemical 
component in a given cell in the simulation box. 
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4 Bertram et al. 


2.4 The A-variance method 


We use the A-variance tool developed by [Ossenkopf et al.| 


( 2008a|p] The A-variance method measures the variance in 


a structure 5'(r) (in our case the maps of centroid veloci¬ 
ties, integrated intensities and column densities) on a given 
spatial scale by filtering the dataset with a spherically 
symmetric up-down-function of size 1. It is given by 


= (('S'(r) * 


( 5 ) 


where the average is computed over all positions r = {x^y) 
on the sky. The symbol * stands for a convolution and 
describes the filter function. In this paper, we use a Mexi¬ 
can hat with a diameter ratio of 1.5. However, we have also 
analyzed the impact of the hlter function and the diame¬ 
ter ratio on the slope values of the spectra, e.g. by using a 
French hat with a diameter ratio of 3.0. We find that our 
results do not significantly depend on the specific choice of 
the filter function and the diameter ratio (see Appendix [b]). 
The differences between the individual filter functions are 
described in more detail in Ossenkopf et al.| ( |2008a| ) . 

In this study we compute the A-variance for each pos¬ 
sible line-of-sight direction x, y and z. As shown by |Esquivel| 
& Lazarian (2005) and Burkhart et al. (2014), the statistics 


of velocity centroids are very sensitive to the direction of the 
magnetic field in the regime of sub-Alfvenic Mach numbers. 
Although we use a weak magnetic held in the ^-direction, 
the turbulence in our simulations is trans-Alfvenic or mildly 
super-Alfvenic (see Table [^. The held lines are essentially 
dragged along with the turbulent how, with the result that 
the turbulence remains approximately isotropic. Hence, we 
do not hnd signihcant variations of the A-variances along the 
different directions. We therefore average all A-variances of 
the three line-of-sights. Finally, power-laws of the form 


CF\{t) OC 


( 6 ) 


were ht to the A-variance spectra, where a denotes the slope 
of the power-law. To calculate the scaling exponents, we 
use a htting range from 1/10 to 1/4 of the total box size 
{0.1 < i/D < 0.25), as constrained by |Federrath et ah (2010) 
and [Konstandin et al. ( |2012| . For a box with D = 512 grid 
cells for each side this translates to 51 and 128 cells, cor¬ 
responding to a physical scale of ~ 2 — 5 pc in the simula¬ 
tion domain. Extending the htting range to scales above this 
limit is complicated since the simulations are driven on large 
scales, which would signihcantly bias our results. Addition¬ 
ally, the scales below our limit are inhuenced by the numeri¬ 
cal resolution and the bottleneck effect ( Kritsuk et al.|20'07 


Konstandin et al.|2015 ), which is an accumulation of kinetic 


energy caused by the viscosity of the huid before it dissipates 
into thermal energy. Eurt her more, as demonstrated by |Kon-| 
standin et al. (2015), we also note that a clear and unique 


identihcation of the htting range in numerical simulations 
with our resolution is complicated. Thus, the specihc choice 
of the htting range introduces some degree of uncertainty in 
the A-variance slopes, which we conservatively estimate to 
be of the order of ± 0 . 1 . 

Eig.[^ shows an example of a A-variance spectrum com¬ 
puted using maps of centroid velocities inferred from our 

^ http://hera.phl.uni-koeln.de/ ossk/Myself/deltavariance.html 


n300 run. The velocities are weighted by the total density 
held and the spectrum is averaged over all available time 
snapshots. We also show the htting slope a for this spec¬ 
trum. This slope can be translated into a linewidth-size rela¬ 
tion, given by cfa{^) oc P^ with a scaling exponent 7 = a/2. 
The latter equation relates the velocity huctuations to their 
spatial scale, which is often referred to as Larson’s law ( |Lar-| 
son||1981 ). This relation suggests that turbulence plays an 
important role in the process of star formation and predicts 
a turbulent energy cascade within the cloud, as proposed by 


Kolmogorov (1941). Hence, the slopes 7 characterize the tur¬ 
bulent velocity hierarchy in our clouds and can be directly 
compared to values derived from spectral observations of 
MCs, typically ranging from 7 0.2 — 0.9 (see, e.g. Larson 


Brunt||2004 

|Heyer et al.| 2006| 

Sun et al.||2006| Hily-Blant 

et al.|20081 

Roman-Duval et al. 

20111 ISchneider et al.|20111 


2014| [Johnston et al.||2014| ). Eor further reading about the 


theory of turbulence in astrophysics and huids, we refer the 


reader to 

Burgers 

(1948), Benzi et al. (1993), Erisch & Kol- 

mogorov (1995|), 

She & Leveque (1994), Dubrulle (1994), 

Boldyrev et al.|(2002) or Mac Low & Klessen (2004). 


Eurthermore, the scaling exponents a of the A-variance 
spectra are related to the power-law exponents /3 of the cor¬ 


responding Eourier spectra via a — (3—1 for 3D data (Stutzki 
et al.|19^ ), where the power-law spectrum in Eourier space 
is defined as 


P{k) oc 


( 7 ) 


In this equation, k — 27r/^ denotes the wavevector. In anal¬ 
ogy, for a power spectrum of a 2 D image, the A-variance is 
related to the power spectrum via a = /3 — 2 . 

We also note that there are many other techniques to 
measure structural density and velocity fluctuations as a 
function of spatial scale, e.g. by computing structure func- 

Eederrath et al.|2010 ) 


tions (see, e.g. Hily-Blant et al. 


2008 


or the spectral correlation function ( Rosolowsky et al.|1999 ), 
by using the velocity channel analysis (Lazarian & Pogosyan 


2000 2004) or the Principal Component Analysis (Heyer & 
Schloerb||1997 Brunt Sz Ileyer|| 20 Q 2 a| Bertram et al.|2(3T4 ). 


3 RESULTS 

We perform numerical simulations and apply radiative 
transfer post-processing to our data in order to analyze the 
influence of chemical inhomogeneities and optical depth ef¬ 
fects on the A-variance analysis. Tablegives an overview of 
our numerical models. Mass- and volume-weighted quanti¬ 
ties are defined via (/)mass = X] f P^V/ P^V and (/)voi = 
Y /dU/ Y ^^5 respectively. Eor more information about the 
H 2 and CO distributions produced in this kind of turbulent 
simulation as well as integrated intensity and column den¬ 


sity PDEs, we refer the reader to |Glover et al.| (pOTOf and 
Shetty et al.| ( 2011 a). 


Eig. 1^ and show the column density and integrated in¬ 
tensity maps as well as the centroid velocity maps computed 
via equations §, § and (0 for all models and chemical 
components. Regarding Eig.j^ we hnd that the total and H 2 
density models show similar structures on all spatial scales, 
which is because most of the hydrogen is in molecular form 
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n30 nlOO n300 



Figure 2. Logarithmic maps of column density as well as velocity-integrated intensity along the 2 ;-direction for our n30 (left column), 
nlOO (middle column) and n300 models (right column). From top to bottom we show the different chemical components: total, H 2 and 
CO column density as well as the integrated intensity of and in the J = 1 —)■ 0 transition. Each side of the simulation domain 

has a length of 20 pc. Note that the velocity field of the n30 model uses a different turbulent random seed than the nlOO and the n300 
model. Furthermore, we caution the reader that our color bars use a different scaling in the individual plots. 
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n30 


nlOO 


n300 
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t»r 


Figure 3. Same as Fig.[^ but with maps of centroid velocities (CV). Black areas in the and intensity map of the n30 model 

denote regions where the brightness temperatures are zero along the line-of-sight. We mask these regions, because no proper centroid 
velocities can be computed there. 
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Model name 

n300 

nlOO 

n30 

Mean density [cm“^] 

300 

100 

30 

Resolution 

5123 

512^ 

5123 

Box size [pc] 

20 

20 

20 

^end [Myr] 

5.7 

5.7 

5.7 

(Ms) 

10.6 

6.8 

5.1 

(Ma) 

1.5 

1.1 

1.0 

(^pKp) 

3.0 

3.0 

4.6 

(^H 2 )mass 

0.98 

0.78 

0.61 

{77H2>vo1 [cm-3] 

140 

37 

9 

(77H2)mass [cm 

1456 

447 

264 

(77co>vol [cm-3] 

2.2 X 10-2 

1.7 X 10-3 

1.4 X 10-4 

{ucO>mass [cm“^] 

0.3 

4.3 X 10-2 

2.3 X 10-2 

(F>voi [K] 

35 

68 

223 

(T)mass [K] 

13 

26 

57 

(iVH2> [cm-2] 

8.6 X 1021 

2.3 X 1021 

5.3 X 1020 

(Nco) [cm 2] 

1.3 X IOI 8 

1.1 X lOi'^ 

8.7 X 1013 


Table 1. Overview of our different models with some character¬ 
istic values measured for the last time snapshot. From top to bot¬ 
tom we list: mean number density, resolution, box size, time of the 
last snapshot, mean sonic Mach number, mean Alfven Mach num¬ 
ber, ratio of density dispersion and mean density, mass-weighted 
mean abundances of H 2 (i.e. the percentage of atomic hydrogen 
that has been converted to H 2 , see [Glover fc Mac Low|[2Qll| >, 
mean volume- and mass-weighted H 2 and CO number densities, 
mean volume- and mass-weighted temperature and mean column 
densities of H 2 and CO. 


at this time (see also Table [^. The CO column densities 
also trace similar structures, but span a much wider range 
of values, demonstrating that carbon monoxide has very low 
abundances along the low column density LoS. Furthermore, 
the and intensity maps also largely reflect the 

distribution of the carbon monoxide gas. We find that the 
intensity maps are much smoother and span a smaller range 
of values (4 — 7 orders of magnitude in integrated intensity 
compared to 7—18 orders of magnitude in column density), 
which is due to the fact that the J = 1 ^ 0 line of 
is easily excited and can be bight even in low-density cloud 
regions. In the case of the intensity maps, the abun¬ 

dance is lower by a factor of i?i 2 /i 3 = 50 compared to the 
intensity models. Nevertheless, is optically thin 

and so the peaks in the maps coincide with those of 

the CO density models. 


3.1 A- variance analysis of the CV maps 

We compute the A-variance for all models (n300, nlOO and 
n30) and all chemical tracers: the total density, H 2 density, 
CO density, ^^CO and ^^CO intensity. We average all spec¬ 
tra from 3 snapshots in time (with 3 line-of-sight directions 
each) where we can assume both the chemistry and the tur¬ 
bulence to be in a stationary and converged state. Fig. 
shows A-variance spectra as a function of spatial scale for 
all density models and chemical tracers, computed for the 


centroid velocity (CV) maps, as introduced in Section 2.3 
Table [^summarizes all slope values ol obtained from a y -fit 
for the A-variance of the velocity field, cr^(^) oc which 
are also listed in each plot in Fig. 

We generally find significant differences between the 
various models and tracers. In our n300 model, we obtain 


Computed on CV maps: 


L[pc] 




L/D 

Figure 4. A-variance spectra as a function of the spatial scale, 
averaged over all available time snapshots and the different lines- 
of-sights y and z. The spatial scale is normalized by the total 
box size. From top to bottom: spectra for our three different den¬ 
sity models, i.e. n300, nlOO and n30, computed for the centroid 
velocity (CV) maps. In each plot we show the A-variance spectra 
for our various chemical components, i.e. for the total density, 
H 2 and CO density as well as for the ^^CO and ^^CO intensity. 
We use a fitting range from 1/10 to 1/4 of the total box size 
(0.1 < i/D < 0.25), as constrained by [Federrath et al.j ( |20fo| 
and [Konstandin et al.j ( |2012| >, which is indicated by a horizontal 
dashed line and is the same in each density model. The different 
power-law functions in the fitting range are indicated with a black 
solid line on each spectrum. Furthermore, we list the slopes a and 
their errors from a each plot. 
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Total density 

a 

H 2 density 

a 

CO density 

a 

intensity 

a 

intensity 

a 

n300 

1.10 ± 0.02 

1.06 ± 0.02 

0.66 ±0.03 

0.75 ±0.03 

0.73 ±0.02 

nlOO 

0.83 ±0.01 

0.70 ±0.01 

0.28 ± 0.01 

0.30 ±0.01 

0.31 ±0.02 

n30 

1.29 ±0.01 

1.27 ±0.03 

0.86 ± 0.02 

0.85 ±0.02 

0.84 ±0.02 


Table 2. Spectral slope of the A-variance spectrum, computed using centroid velocities weighted by the indicated quantity (see also 
Fig. The errors are computed by the fitting method. The slope a is related to the A-variance via cr^(^) oc in a linear regime. 
The values a can also be used to compute spectral slopes 7 = a/2 for a linewidth-size relation (7/\{t) oc readily comparable to 
observational measurements. 


similar slopes for the total density and the H 2 density (see 
Table [^. Regarding the spectra of those two cases, we find 
an excellent agreement within the values of cr^ for both mod¬ 
els, as shown in the top plot of Fig. For the low-density 
model n30 (bottom plot of Fig. [^, we obtain a large dis¬ 
crepancy between the values of total density and the H 2 
density model. This is because in this simulation the frac¬ 
tion of molecular gas is much smaller than in the higher 
density runs and so H 2 no longer follows the total gas den¬ 
sity ( Glover et al.||2010 ). The H 2 density is therefore more 
inhomogeneous than the gas density, and as a result the cen¬ 
troid velocities weighted by H 2 fluctuate more on all scales 
than those weighted by the total density. The nlOO model 
makes up an intermediate case between the n300 and the 
n30 density models. In this model, the correlation between 
the total and the H 2 density is worse than in the n300 case, 
but still better than in the n30 model. 

Comparing the total and H 2 density models to the var¬ 
ious CO tracers, we find that the slopes of the former are 
significantly steeper compared to the slopes of the latter, 
independent of the density (see Table [^. However, the A- 
variance spectra of the various CO tracers show a good 
agreement with each other over nearly all spatial scales. Fur¬ 
thermore, we always find the a\ values of the CO tracers to 
be significantly larger than those of the total and H 2 density. 
This is because CO is mainly located in denser regions of the 
cloud (see, e.g. Shetty et al. 201 la,b), leading to higher den¬ 
sity contrasts compared to the total and H 2 density and thus 
to larger variances in velocity space. 


3.2 A-variance analysis of the intensity and 
column density maps 

Beyond the A-variance analysis of centroid velocities, we can 
also apply the same method to any other quantities defined 
on the x-y plane. In particular, it is possible to carry out a 
similar analysis for the total, H 2 and CO column densities 
and the ^^CO and ^^CO integrated intensities. The results 
of this analysis are illustrated in Fig. In contrast to Fig. 
1^ we normalize each curve by the corresponding A-variance 
value ctq measured at an arbitrary spatial scale of 5% of the 
total box, in order to better compare the spectra computed 
for the intensity and column density maps with each other. 
Furthermore, we also show a horizontal dashed line at a 
value of cr^/cTo = 1 in each plot. Table summarizes all 
slope values a obtained from a for the individual A- 

variance spectra. 

As for the CV maps, we find significant differences be¬ 
tween the spectra of the various models and tracers com¬ 


puted using the maps of integrated intensity and column 
density directly. In the high-density n300 model, we find 
a good correlation between the total density and the H 2 
density (see also the slopes in Table , for the same rea¬ 
sons as described in Section [3.1| Similarly, the low-density 
n30 model shows again the largest discrepancy. This is be¬ 
cause the fraction of molecular gas is significantly smaller 
than in the n300 simulation, leading to different values of a. 
The nlOO model again defines an intermediate case between 
the n300 and the n30 density runs. In this model, the H 2 
density better correlates with the total density than in the 
low-density n30 run, but still worse than in the high-density 
n300 run. 

Regarding the simulations of the various CO tracers, we 
obtain very different scaling properties compared to the to¬ 
tal and the H 2 density simulations, depending on the mean 
density. In the high-density n300 model, we find a good cor¬ 
relation between the two CO intensity cases, while the slopes 
of the total and CO density models are significantly steeper. 
However, both the total and H 2 density as well as the var¬ 
ious CO tracer models show the same trend in over all 
spatial scales, having a > 0. This situation changes in the 
low-density n30 model, and we find strong differences be¬ 
tween the CO tracers and the total and H 2 density runs over 
all spatial scales. While the structure in the computational 
domain on larger scales for the total and H 2 column density 
can still be described by a slope a > 0, the A-variance spec¬ 
tra for the CO tracers have negative slopes in this density 
run (see Table . A similar situation as in the low-density 
n30 simulation can be seen in the intermediate nlOO density 
model. While we again measure a > 0 for the total and H 2 
density cases, we obtain a < 0 for the CO tracers, meaning 
that the spatial distribution of observable carbon monoxide 
is completely different compared to the distribution of the 
H 2 gas in the n30 and the nlOO density simulations. 

Looking at Fig.|^ we see that the reason that we obtain 
negative values for a for the CO tracers in the low-density 
runs is that the corresponding A-variance spectra turn over 
at relatively low spatial scales. This indicates that the CO in 
these simulations has a characteristic spatial scale of approx¬ 
imately 2% of the total box size, corresponding to around 
0.4 pc. This is consistent with what we see in the CO column 
density projections and integrated intensity maps shown in 
Fig. in In the low-density n30 simulation, CO is sufficiently 
well shielded to resist photodissociation only in the dense 
filaments formed by turbulent compressions. These struc¬ 
tures are thin and contribute strongly to the A-variance on 
scales comparable to their width. Moreover, the fact that 
this characteristic scale is roughly the same in all CO spec- 
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Total density 

OL 

H2 density 

OL 

CO density 

OL 

intensity 

OL 

^^CO intensity 

OL 

n300 

0.34 ±0.02 

0.32 ±0.02 

0.37 ±0.02 

0.23 ± 0.02 

0.22 ±0.01 

nlOO 

0.20 ±0.01 

0.10 ±0.01 

-0.18 ±0.01 

-0.04 ±0.01 

-0.06 ±0.01 

n30 

0.51 ±0.01 

0.39 ±0.03 

-0.98 ± 0.04 

-0.35 ±0.04 

-0.72 ±0.01 


Table 3. Same as Table but for the A-variance computed directly from maps of integrated intensity and column density, as defined 
in equation |3) and lUl. The corresponding A-variance spectra are shown in Fig.js] 


tra demonstrates that the radiative transfer post-processing 
does not significantly change the characteristic spatial cloud 
signatures in the A-variance analysis. Furthermore, is 

optically thick and thus we always find slightly steeper a 
slopes for the intensity tracers compared to the opti¬ 
cally thin intensity tracers. 

3.3 Variation of the abundance of 

We also analyze the variation of the abundance of ^^CO, 
because our simulations only follow the chemistry of ^^CO, 
but not of its isotope ^^CO. Therefore, we produce a set of 
emission maps, using a spatially varying i?i 2 /i 3 gen¬ 
erated following the prescription given in Section |2.2| and 


Szucs et al. (2014). After applying the radiative transfer 


post-processing to these maps, we compute the maps of the 
CV and the integrated intensity and compare the A-variance 
spectra with each other. The results agree with those as de¬ 


scribed in Bertram et al.| (2015a) for the centroid velocity 
increment structure functions. Although we encounter slight 
variations on small scales for both the CV and the inte¬ 
grated intensity maps, the spectra of the A-variance do not 
change signihcantly as we vary i?i 2 /i 3 - Thus, the use of a 
constant ^^CO/^^CO ratio is sufficient for obtaining proper 
A-variance spectra for ^^CO. 


4 DISCUSSION 

4.1 Interpreting the A-variance spectra computed 
using the CV maps 

Regarding the spectra computed on the CV maps in Fig. 
1 ^ we hnd a good correlation between the total density and 
the H 2 density in the n300 model. This is because ~ 98% 
of the initial atomic hydrogen is in molecular form at this 
time. Thus, the spectra of the A-variance for the CV maps 
are primarily dominated by the H 2 mass in this case. Vice 
versa, if we analyze the CV spectra of the n30 model, we hnd 
a larger discrepancy between the total density and the H 2 
density, which is because only ~ 61% of the initial atomic 
hydrogen is in molecular form at this time. The nlOO model 
is an intermediate case, with this value being ~ 78%. Fur¬ 
thermore, we obtain a good correlation between the various 
CO tracer components. This indicates, that the turbulence 
statistics are similar for ah three cases and that the impact 
of the radiative transfer post-processing on the A-variance 
analysis is limited. This is similar to the conclusion pre¬ 


sented in Bertram et al. (2015a) based on the analysis of 


for the different CO tracer models (see the values in Ta¬ 
ble]^. This indicates that these components have a higher 
relative amount of structures on larger scales compared to 
the different CO tracers. Furthermore, as shown by |Glov^ 
et al. (2010) and She tty et al.| (2011a,b), CO is primar¬ 
ily a good tracer of dense and very compact regions in a 
cloud. Thus, this leads to less turbulent velocity structures 
on larger scales and hence to hatter slopes compared to those 
of the total and H 2 density. However, we hnd that the dif¬ 
ferent slopes a of the various CO tracers underestimate the 
slopes of the total and H 2 density by a factor of ~ 1.5 — 3.0 
(see Table [^. 

Furthermore, we caution the reader that it is difficult 
to infer a clear dependence of the slopes on the mean ISM 
density. It is likely that the statistical measures we derive 
from our numerical simulations are also sensitive to the spe- 
cihc realization of the turbulent velocity held. Since we are 
studying hows which are driven on large scales, variance ef¬ 
fects can become important, and the statistical properties 
depend on the random orientation of the turbulent modes 
as well (see, e.g. Klessen|2Q00 Klessen et al.|2000 ). Hence, in 
order to obtain slopes that properly converge with the den¬ 
sity, we speculate that a large number of simulations with 
various turbulent realizations would be needed in order to 
calculate reliable average values (see, e.g. the PCA analysis 
of the statistics using different realizations of the turbulent 
velocity held in Bertram et al.|[2014| ). This is prohibited by 
the high computational costs of the individual simulations, 
and hence we only focus on one example in this paper, which 
is enough to illustrate basic trends of the A-variance statis¬ 
tics. 

We also note that a direct comparison of our 2D A- 
variance statistics to 3D turbulence models is difficult, since 
ah our 2D maps are a complex convolution of the 3D density 
held (or the brightness temperatures) with the 3D veloc¬ 
ity held. Various different physical processes can inhuence 
the CV statistics, as shown by previous studies. For exam¬ 
ple, Lazarian V Pogosyan (2004), Burkhart et al.| (2013a) 


and Burkhart et al. (2014) analyzed the impact of the sonic 
Mach number on the CV statistics, hnding that it can sig- 


nihcantly alter the results. Moreover, |Lazarian & Esquivel 

(2003), Ossenkopf et al. 

(2006), Esquivel et al. 

(|2007|),|Hily- 

Blant et al. (2008) and 

Eederrath et al. (2010 

) studied the 


CV increment structure functions. 

In general, we hnd that the slopes a for the total and 
H 2 density models are signihcantly steeper than the slopes 


effects of the turbulent driving as well as temperature and 
density huctuations on the CV statistics, also hnding sig- 
nihcant differences in the statistics and in the inferred CV 
slopes. Nevertheless, although a direct comparison of the 2D 
statistical quantities to 3D measures is complicated, we can 
safely use the 2D A-variance analysis in order to work out 
statistical trends measured in the spectra of the individual 
models. 
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CO intensity, a = -0.04 V- 0.01 
CO intensity, a = -0.06 V- 0.01 



L/D 

Figure 5. Same as Fig.^ but computed using maps of integrated 
intensity and column density. In this plot, we normalize each curve 
by the corresponding value ctq measured at an arbitrary spatial 
scale of 5% of the total box, in order to better compare the spectra 
computed for the intensity and column density maps with each 
other. For a better visualization, we also show a horizontal dashed 
line at a value of ct^/ctq = 1 in each plot. Furthermore, for clarity, 
we do not plot the error bars of the individual spectra here. The 
slopes a and their errors from a listed in each plot for 

all our chemical components. 


4.2 Interpreting the A-variance spectra computed 
on the intensity and column density maps 


In Section |3.1| and |3.2| we already justified the correlation 
between the total and the H 2 density runs in all our density 
models and compared the different slopes a from Table 
with each other. Furthermore, we established that the ob¬ 
served CO gas distribution in the low-density n30 model 
does not reflect the spatial distribution of molecular hy¬ 
drogen very well. This is because CO is mainly located in 
the dense filaments, while H 2 is more space-filling and dis¬ 
tributed over the total cloud, owing to its greater ability to 
resist photodissociation. Consequently, the A-variance spec¬ 
tra of the CO tracers peak at a scale corresponding to the 
width of these structures. Hence, if we apply the A-variance 
analysis to the maps of integrated intensity or column den¬ 
sity, we can obtain important information about the turbu- 
lently created high-density peaks within the MC (see also 
the discussion in Ossenkopf et al.||2001 ). However, we also 
note that the physical connection of our slopes a from Table 
to observational measurements is complicated, since ob¬ 
servations typically probe smaller spatial scales, which our 
simulations are not sensitive to due to the limited numerical 
resolution. 


If we compare the CO tracers with the total and H 2 den¬ 
sity in the n300 model, we find that the carbon monoxide 
better reproduces the statistical trends of molecular hydro¬ 
gen in this high-density run compared to the low-density 
n30 model. In the n300 run, we generally obtain a > 0 for 
all chemical components, i.e. we find more cloud structures 
on larger spatial scales. In contrast to the low-density n30 
model, this is because carbon monoxide is not only confined 
to small dense filaments, but instead is distributed over the 
whole molecular cloud. In particular, a significant amount 
of CO gas can also be observed between the numerous dense 
cloud regions. Consequently, we see the A-variance spectra 
peaking at the largest scales, consistent with observational 


efforts and previous work on this field (see, e.g.|Stutzki et ah 


1998 

Ossenkopf et al. 

[2001 

IBensch et al. |2001| 

Sun et al. 

2006 

Schneider et al. 

poll 

Russeil et al.||2013 

|Alves de 


Oliveira et al.||2014[ [Elia et al.||2014|). 


However, regarding the general trend of the slopes a 
from the n30 to the n300 via the nlOO model, we find that 
the values of a for the CO tracers change sign above a critical 
density of ~ 100 cm“^. This is about the number density at 
which all values a become positive. Thus, we conclude that 
carbon monoxide traces our total cloud structure well only 
if the average cloud density lies significantly above a critical 
threshold of ~ 100 cm“^. If the mean density in the cloud is 
significantly smaller than this limit, the observable CO gas 
does not properly trace the statistical properties of the H 2 
gas in the cloud. We speculate that one may also see a similar 
switch from CO tracing all of the structures to only the 
dense cores and filaments if one increases the incident field 
strength or decreases the metallicity, as in both cases, this 
makes it much easier to photodissociate the diffuse CO (see, 
e.g. Glover & Clark 2012; Bertram et al. 2015b, submitted). 


On the contrary, if the mean density in the cloud is too 
high, CO becomes optically thick and so we suspect that 
there should also exist an upper density limit above which 
CO does not properly trace structures of the cloud anymore. 
We leave such an analysis for further investigations. 
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4.3 Comparison of the A-variance to other 

statistical tools applied to our simulation data 

We have previously applied several other statistical methods 
to the same set of simulation data in order to study the 
structural behavior of MCs. For example, [Bertram et al.| 
(2014) applied Principal Component Analysis (PCA) to the 
same data used in this paper, while Bertram et al.| (2015a) 
carried out a similar study using centroid velocity increment 
structure functions (CVISF) and Fourier spectra. 

Comparing the results of the different methods with 
each other, we generally find consistent results between 
the different statistical analysis methods. For example, the 
slopes of the CVISF for the total density and H 2 density 
models are also significantly steeper compared to the vari¬ 
ous CO tracers, which is in good agreement to the results 
of the A-variance analysis in this paper. The same holds for 
the relative scaling of the different CO models, which is the 
same in the analysis of the CVISF and the A-variance. In¬ 
terestingly, we also find a close relation between the PCA 
structure analysis and the spectra of the A-variance for the 
CO tracers. The PCA method does not find any structures 
on larger scales for the CO tracers in the n30 model as well 
as for some single CO tracers in the nlOO model. This re¬ 
sult can also be reproduced in the different CO spectra of 
the A-variance analysis, i.e. where the gradient of the indi¬ 
vidual spectra become negative at a characteristic scale of 
~ 0.4 pc, as described in Sections |3.2| This is the situation 
where small clumpy structures of CO gas dominate the over¬ 
all composition of carbon monoxide in the MC. In this case, 
the A-variance shows that the most dominant CO structures 
in the simulation domain can be found at small scales, lead¬ 
ing to completely missing principal components on scales 
larger than ~ I pc, as presented in Fig. 3 in [Bertram et al.| 
(2014). Hence, applying the A-variance analysis to maps of 


of a cloud midway between the diffuse and the molecu- 


integrated intensities or column densities gives an idea about 
characteristic spatial scales in the cloud of interest. The A- 
variance thus provides a good statistical tool in order to 
study the relative gas distribution on various cloud scales. 

These results indicate that all these various statistical 
methods are connected to each other and that characteristic 
structural properties of the MC should be traced by each 
of them individually. Hence, all various statistical meth¬ 
ods have proven to yield reasonable results for the struc¬ 
ture analysis of MCs. Nevertheless, the advantage of the 
A-variance method is that it is fast and easy to implement, 
while the computation of structure functions and principal 
components is more expensive. Thus, the A-variance anal¬ 
ysis provides a useful and adequate tool in order to quickly 
study the internal structure of a cloud. However, the advan¬ 
tage of the structure function analysis is that their results 
can be easily compared to theoretical models of turbulence. 


4.4 Previous studies about A-variances computed 
using CV maps 

Several studies tried to compute the A-variance spectra of 
CV maps in the past and estimated turbulent slope values a 
from observa tional measurements. For example, [Ossenko']^ 


et al. 


(2008b) used optically thick ^"^CO (J = I ^ 0) maps 


of the Polaris Flare to compute centroid velocity maps and 
to estimate the slope a. The Polaris flare is an archetype 


lar phases ([Heithausen & Thaddeus 1990 Meyerdierks & 


Heithaus^|I996| [Falgarone et al.[[I998| [Andre et al.[[20I0[ 

Miville-Deschenes et al.[[20I0[ ). It is supposed to have a low 

average density and thus can be compared to our n30 low- 
density model. Adopting the same ^^CO tracer, we And a 
value of a — 0.85 ± 0.02, in agreement with the estimate of 


OL 0.9 found by Ossenkopf et al. (2008b). 


Moreover, Ossenkopf & Mac Low (2002) computed the 


slope OL for hydrodynamic supersonic simulations driven 
at different wavenumbers. They And one power-law range 
for all models and obtain ol ~ I.O, which fits into our 
range of power-law slopes ol from 0.8 — 1.3 for the total 


density models. Furthermore, Federrath et al. (2009) com¬ 


puted A-variance slopes for numerical simulations with both 
solenoidal and compressive forcing, finding a similar range 
of OL values from 0.8 — 1.4 for the turbulent velocity held. 
However, these simulations only use supersonic isothermal 
turbulence, while our runs include more complex physics, 
e.g. a chemical network, heating and cooling, various initial 
number densities or the coupling to the radiation held. Thus, 
if we compare the total variations of slopes a in their and in 
our models, we And that the influence of the different forcing 
methods on the slopes ol in their simulations is large, while 
the impact of our additional physical effects (varying den¬ 
sity, optical depth effects, etc.) on the slope values remains 
comparatively small. 

4.5 Previous studies about A-variances computed 
using intensity and column density maps 

As well as computing A-variance spectra for maps of cen¬ 
troid velocities, several studies also applied the A-variance 
to maps of integrated intensity or column density. For ex¬ 


ample, Bensch et al. (2001) computed A-variance spectra 
on maps of velocity-integrated intensity for and 

(J — 1 _ ^ 0) for various MCs in the Galaxy. Stutzki 


et al. (1998) analyzed the A-variance of an observed 


(J = I ^ 0) image of the Polaris flare as a whole, while 


Ossenkopf et al. (1998) studied the intensity map of one of 


its subclouds, MOLD 123.5+24.9. 

Comparing all the various A-variance spectra in [Ben-[ 


sch et al. (2001) or Stutzki et al. (1998) with each other. 


we And a value a > 0 for each of them in the given fitting 
range, i.e. the spectra increase with increasing spatial lag. 
At first sight, this is in contradiction to the results that we 
obtain in this study, where we measure a < 0 for the vari¬ 
ous CO tracers in the n30 and nlOO model, given our fitting 
range. However, this is primarily due to the fact that our A- 
variance spectra peak at roughly the scale of small carbon 


monoxide clumps in the low-density clouds (see Section 4.2). 
Accordingly, we And more cloud structures on the scales of 
those localized CO structures, leading to negative ol slopes 
in Figure A similar effect can be seen in the A-variance 
analysis of ^^CO (J = I ^ 0) maps of the outer Galaxy 
shown in Figure 3 of [Stutzki et al. (1998). On small scales, 
a positive power-law slope is recovered, but on large scales 
a becomes negative as one starts to smooth on scales larger 
than the typical size of the observed MCs. The same situa¬ 
tions aFoillustrated in the A-variance spectra in [Schneider[ 


et al. (2011) for the ^^CO (J = I ^ 0) molecular line survey 
of Cygnus X as well as in a ^^CO (J = 2 ^ I) study of the 
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Perseus cloud in Sun et al. (2006), which peak at intermedi¬ 
ate cloud scales. However, a A-variance analysis of the small 
sub-regions NGC1333 and L1455 in the Perseus molecular 
cloud complex shows that a is always positive (see Fig. 6 
Sun et ar]|2006 ), which is in agreement with our finding 


that a positive slope a might be related to a A-variance 
analysis at lower spatial cloud scales, probing only very lo¬ 
calized CO structures. Moreover, [Alves de Oliveira et al.| 
(2014) analysed the Chamaeleon molecular cloud complex 
using data from the Herschel Gould Belt Survey. This com¬ 
plex encompasses three MCs with different star formation 
histories. One of them is a quiescent cloud, which should 
be best comparable to our numerical simulations, since we 
neglect the effects of self-gravity. The A-variance analysis 
of this MG also shows a clear break at mid to large scales 
(see cloud Cha III in Fig. 5 in Alves de Oliveira et al.|20T4 ), 
also displaying a characteristic spatial scale in the column 
density structure. 

In our high-density n300 run, the CO gas is distributed 
over the whole MG and not only confined to small dense frag¬ 
ments. In this case, we also find cloud structures on larger 
spatial scales, leading to positive a slopes in Figure 


4.6 Model limitations 

Since we are running numerical simulations, we have to keep 
in mind that our runs are subject to various physical simpli¬ 
fications. In the current analysis, we focus on the impact of 
the chemistry and the opacity of the gas on our results, thus 
neglecting other physical processes that could bias our anal¬ 
ysis. In first instance, we do not account for self-gravity and 
thus also do not model star formation or stellar feedback. 
More specifically we do not account for stellar radiation, SN 
feedback and other physical processes. We also do not in¬ 
clude any large-scale dynamics, e.g. spiral arms or galactic 
rotation, although these are unlikely to be important on the 
20 pc scale studied here. Nevertheless, we can infer useful 
information about how the chemical composition of the gas 
and the opacity affect the A-variance analysis. We also note 
that our results depend only weakly on the resolution and 
that this concerns mostly GO as a tracer molecule, as we 
show in Appendix For future investigations, we want to 
analyze simulations that span a wider range of physical pa¬ 
rameters, e.g. with different levels of the external radiation 
field, varying metallicities or additional physical processes, 
in order to find out how they affect the statistics. 


• We compute A-variance spectra of maps of centroid ve¬ 
locities and fit a power-law cr^(^) oc in order to charac¬ 
terize the properties of the turbulent hierarchy in the MGs. 
This power-law can be translated into a linewidth-size rela¬ 
tion, i.e. crA(^) oc P with 7 = q;/ 2 , readily comparable to 
slopes that can be derived from spectral observations. We 
find the slopes a of both the total and H 2 density models 
to be significantly steeper than the slopes of the different 
GO tracers, which underestimate the former by a factor of 
~ 1.5 — 3.0 (see Section 3.1). 

• The slopes a derived from the GV maps range from 0.8 
to 1.3 for the total and H 2 density, while a for the various GO 
tracers range from 0.3 to 0.8 (see Section 3.1). However, we 
also note that the specific choice of the fitting range might 
cause further variations of the slopes by ±0.1. 

• Although we find slight variations between the different 
slopes a for our various GO models obtained from the GV 
maps, the impact of the optical depth effects on the spectra 
computed on maps of centroid velocities remains limited (see 
Table 1 ^. 

• This is different in the case of the integrated intensity 
and column density. The A-variance computed from these 
maps is strongly affected by optical depth effects. The GO 
tracers exhibit a very different spatial scaling behavior com¬ 
pared to the total and H 2 density models (see Table and 
Figure]^. 

• We report a critical number density threshold of ~ 
100 cm“^ at which the spectral slopes a of the GO trac¬ 
ers switch sign for the A-variance of integrated intensity 
and column density maps. We conclude that carbon monox¬ 
ide traces the total cloud structure well only if the average 
cloud density lies significantly above this threshold. If the 
mean density in the cloud is significantly smaller than this 
limit, the observable GO gas does not properly trace the 
statistical properties of the H 2 gas in the cloud (see Section 


3.2). 


• The A-variance spectra computed on maps of inte¬ 
grated intensity and column density provide a useful statis¬ 
tical measure in order to infer important information about 
the distribution of gas within a cloud. We also argue that 
peaks in the A-variance spectra correspond to characteris¬ 
tic scales of the morphological structure of the system (see 
Section 3.2). 

• Our findings are consistent with previous A-variance 
studies using GO line observations or measurements of the 
continuum (see Section 


4.4 


and 4.5). 


5 SUMMARY AND CONCLUSIONS 

We analyzed A-variance spectra of MCs with time- 
dependent chemistry and radiative transfer post-processing 
for models of different initial number densities and chem¬ 
ical components: the total number density, H 2 and CO 
density (each without radiative transfer) as well as ^^CO 
(J = I ^ 0) and ^^CO (J = I ^ 0) intensity (both with 
radiative transfer). In each case, we computed A-variance 
spectra for maps of centroid velocities (GV), integrated in¬ 
tensities and column densities and analyzed the structural 
behavior of MCs in numerical simulations. We report the 
following findings: 
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Figure Al. Same as Fig.^ but with runs of different resolutions 
with 512^ and 256^ grid cells with a fixed initial number density 
of no = 100 cm“^, computed for the centroid velocity (CV) maps. 
We find similar a values for the total density and the H2 density 
models, while the values a differ by up to ~ 20 — 30% between 
the different resolution models for the various CO tracers. 
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Texas Advanced Computing Center, using time allocated as 
part of Teragrid project TG-MCA99S024. Additional sim¬ 
ulations were performed on the kolob cluster at the Uni¬ 
versity of Heidelberg, which is funded in part by the DFG 
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APPENDIX A: RESOLUTION STUDY 

We study the influence of the numerical resolution on the re¬ 
sults of the A-variance. Therefore, we have performed runs 
with resolutions of 512^ and 256^ grid cells and evaluate the 
spectra and slopes for the CV maps for all chemical compo- 


L[pc] 



L/D 

Figure Bl. Same as Fig. but with the A-variance spectra for 
the H2 and CO density models for a fixed initial number density 
of no = 100 cm“^ and a resolution of 512^ grid cells, computed for 
the centroid velocity (CV) maps. In order to analyze the impact 
of the filter function and its diameter ratio on our results, we 
compute the H2 and CO density spectra using both a Mexican 
hat with a diameter ratio of 1.5 and a French hat with a diameter 
ratio of 3.0. We do not find any significant differences in the slopes 
a within the fitting errors if we use a another filter function for 
the A-variance analysis combined with a different diameter ratio. 


nents for a fixed initial number density of no = 100 cm“^. 
The results and their interpretation is the same for all other 
density models. Fig. |A1| shows spectra of the A-variance 
with the corresponding slopes a, in analogy to Fig. The 
fitting range for the 256^ model is downscaled by a factor 
of 2 compared to the 512^ model, i.e. we fit from 25 to 64 
cells in the spatial domain. We find similar a values for the 
total density and the H 2 density models. However, the slope 
values a differ by up to ~ 20 — 30% between the differ¬ 
ent resolution models for the various CO tracers. This is 
because CO is mainly located in dense regions of the cloud 
( Bertram et al.|2014 ), which can be resolved more accurately 
at a higher resolution, leading to significant differences be¬ 
tween the two resolution models. These results agree with 
the results in the resolution study of the structure function 
analysis presented in Bertram et al. (2015a). However, the 
variations caused by the effect of resolution are rather small, 
since a also strongly depends on the specific choice of the fit¬ 
ting range, which might also cause slope variations by about 
±0.1. Nevertheless, we find a similar relative scaling behav¬ 
ior in the two resolution models between the spectra of the 
different chemical components. 


APPENDIX B: COMPARISON OF SPECTRA 
WITH DIFFERENT FILTER FUNCTIONS 

We also study the influence of the filter function and the 
choice of the specific diameter ratio on our A-variance spec¬ 
tra. Therefore, as an example. Fig. |B1| shows the A-variance 
spectra for the H 2 and CO density models for a fixed ini¬ 
tial number density of no = 100 cm“^ and a resolution of 
512^ grid cells, computed on maps of centroid velocities. 
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For each model, we evaluate the spectra using different fil¬ 
ter functions and diameter ratios. In particular, we compute 
the A-variance spectra using a Mexican hat with a diameter 
ratio of 1.5 as well as a French hat with a diameter ratio of 
3.0. In analogy to Section [ffT] we fit a power-law within a 
given fitting range to the spectra and compare the slopes 
with each other, which are shown in Fig. |B1[ Thereby, we 
do not find any significant differences between the slopes de¬ 
rived from spectra with various filter functions and diameter 
ratios within the fitting errors for our models. The individ¬ 
ual normalizations of the spectra are caused by the variable 
choice of the diameter ratio, affecting the computation of 
the A-variance analysis (see, e.g. Ossenkopf et al. 2008a). 
However, the shape of the individual spectra for one chem¬ 
ical model over various spatial scales is also approximately 
conserved. 
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